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Abstract. The use of type la supernovae as distance indicators for cosmology has initiated a search for theoretical arguments 
supporting the empirical calibration methods applied. To this end, as a first step, a sound understanding of the origin of the 
observed diversity in type la supernova properties is needed. Here we present a first systematic study of effects resulting 
from changing some physical parameters of three-dimensional deflagration models of thermonuclear supernovae. In our study 
we vary the progenitor's carbon-to-oxygen ratio and its central density prior to ignition because both properties are not well 
determined by stellar evolution theory and they may change from supernova to supernova. Next we compute for these explosion 
models the nucleosynthesis yields in a post-processing step. This, in addition, allows us to study variations in the progenitor's 
metallicity by means of different 22 Ne mass fractions in the initial composition. We find that the progenitor's carbon-to-oxygen 
ratio and its central density affect the energy release of the models and thus the expansion velocity of the supernova. Moreover, 
we find that changing the metallicity and the central density changes the production of radioactive 56 Ni and thus affects the 
luminosity. In contrast, the carbon-to-oxygen ratio has little effect on the 56 Ni production. Implications of the found variations 
of the explosion energy and the produced 56 Ni mass for the type la supernova diversity are discussed. 
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1. Introduction 

Type la supernovae (SNe la) have become a major tool to de- 
termine cosmological parameters. As a consequence of their 
uniformity of properties and their enormous brightness they are 
suitable to be applied in cosmological distance measurements. 
However, they cannot be claimed to be perfect "standard can- 
dles", as they show a significant intrinsic scatter in their peak 
luminosities as well as other characteristics. Therefore their 
cosmological application rests on empirical corrections of the 
peak luminos i ties ba sed on correlations with other observables 
(e.g. IPhillipsL I1993I) . Only such empirical corrections facili- 
tated distance measurements of SNe la at high redshifts which 
have led to the spectacular conclusion that the expans ion of uni- 
verse currently accelerates jRiess et al.lll998tlPerlmutter et all 
1999). One possibility to technically take this result into ac- 
count is by a cosmological constant in the Einstein equations 
eventually indicati ng a dark energy component of the universe 
(for a review see Leibundgut 2001 ). 

This underscores the need for a theoretical reasoning of 
the correlation of characteristics that are yet established only 
empirically. A theoretical understanding will help to answer 
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questions, such as a possible affliction of calibration procedures 
with evolutionary effects. 

In the astroph y sical standard model (see 
iHillebrandt & Niemevetl l200fi . SNe la are associated 
with thermonuclear explosions of carbon/oxygen white dwarf 
(WD) stars. The optical event is powered by the decay of 
radioactive species (e.g. 56 Ni) produced in the thermonu- 
clear burning. Numerical simulations on the basis of this 
scenario provide an approach to the understanding of cal- 
ibration methods. Recently, there has been much progress 
in the three-dimension al mod el ing of the exp l osion pro- 
cess faillebrandt et all l200fi iReinecke et all l2002alclbt 
iGamezo et all l2003h and the question arises whether it is 
possible to reproduce the SN la diversity by varying the 
initial parameters of such models. This will be addressed in 
the present study where we restrict the survey to so-called 
deflagration models of thermonuclear supernovae which can 
be summarized as follows. 

After ignition near the center of the WD the flame propa- 
gates outward in the subsonic deflagration mode, i.e. it is medi- 
ated by thermal conduction of the degenerate electron gas. This 
outward burning produces an inverse density stratification in 
the gravitational field of the WD star with dense fuel on top of 
hot and light ashes. Consequently, due to buoyancy (Rayleigh- 



2 



F. K. Ropke et al.: Type la supernova diversity in three-dimensional models 



Taylor) instabilities burning bubbles form that rise into the fuel 
leading to shear flows. Kelvin-Helmholtz instabilities generate 
strong turbulence given the fact that the Reynolds number typ- 
ical for the situation is of the order of 10 14 . Resulting from 
this, turbulent eddies decay to smaller scales thereby forming 
a turbulent energy cascade. The interaction of the thermonu- 
clear flame with these turbulent motions is the key feature of 
the deflagration model for SNe la. A laminar flame would burn 
far too slowly to release sufficient energy for an explosion of 
the star. However, the wrinkling of the flame due to turbulence 
and the accompanying flame surface enhancement increase the 
net burning rate and accelerate the flame. This defines the de- 
flagration model of thermonuclea r supernova explosions as a 
problem of turbulent combustion. iReinecke et all d2002bl) and 
IGamezo et alJ d2003l) could show that this model indeed leads 
to an explosion. Whether it reproduces al l aspects of observe d 
supernovae is still no t ful ly explored (e . g.lKozma et all Eo05). 
IGamezo et ail J2004 and Hoflic h et al] d 19981) claim"that a hy- 
pothetical transition from the deflagration mode of flame prop- 
agation to a supersonic detonation needs to be invoked at later 
phases of the explosion. We set aside suc h a transition be - 
cause its physical origin in not understood dNiemeverl fl999>. 
Moreover, even in such a case the initial deflagration stage will 
be essential for understanding the SN la diversity since large 
fractions of the energy and the radioactive 56 Ni (which powers 
the lightcurve) are produced here and nonlinear effects in flame 
propagation are extremely sensitive to the initial conditions. 

The crucial role played by three-dimensional effects in de- 
flagration SN la models calls for multi-dimensional simula- 
tions to study the diversity of such events. Most previous at- 
tempts to unveil the origin of the SN la dive rsity were, however, 
based on one-dimensional mode ls (Bravo etal., 19931 199; 



We restrict this first systematic study to variations of the 
central density, the initial carbon-to-oxygen (C/O) ratio and 
the metallicity of the progenitor just prior to ignition. Our in- 
tention is to test the parameters independently, setting aside 
any realistic evolution of the progenitor syste m. For det ailed 



Hoflic h et allll998tlUmeda et al.lil999illwamoto et allll99 
Domfnguez & Hoflichl l2000t iDomfnguez et all l2000t l200ll) 
These are hampered by introducing free parameters due to in- 
complete description of the relevant physics in addition to the 
initial parameters they intend to study. The description of the 
turbulent mixing as well as the effective flame velocity is not 
inherent in one-dimensional models but rather accomplished in 
a parametrized way. Due to the free parameters empirical one- 
dimensional models are not sufficiently predictive to nail down 
explanations for the diversity of SNe la, but they can never- 
theless provide valuable clues for possible trends. Systematic 
studies based on three-dimensional models overcome the am- 
biguity of the turbulent flame velocity and mixing. By cor- 
rectly modeling these effects, multi-dimensional deflagration 
models contain no tunable parameters and possess a high pre- 
dictive power. However, due to the challenging computational 
demands of three-dimensional models, the available studies of 
initial parameters are very incomplete. 

Applying a simplified setup we present the first systematic 
survey of the impact of initial parameters on three-dimensional 
SN la models. The price of the simplicity (and possibly incom- 
pleteness) of our models is that we cannot set the absolute scale 
of the effects in the presented parameter study. Nevertheless, 
we are able to point out the trends of effects from varying the 
initial parameters. 



progen itor evolution studies we refer to e.g. |N omoto & Iben 
( 198 51) iHernanz et al.ldl988l).lBravo et alJdl99d).IUmeda et al 



Jl999bl) . lLanger et alJ ( Ebool) . and iDommguez et alJ J200ll) . 
Important parameters that are not addressed in this study are 
for instance rotation and the way of flame ignition (see e.g. 
Wooslev etal., 2004). Some effects of the ignition conditions 
on SN la models and nucleosyn thesis have been recently dis- 
cussed bv lTravaglio et al.l (2004). 

In Sect. |2 we describe the numerical schemes we apply to 
model SNe la explosions and the nucleosynthesis, followed by 
a discussion of the parameter space to be explored in Sect. 
The features of the explosion models will be compared in 
Sect. |5] and Sect.[6]describes the nucleosynthetic yields of these 
models. Conclusions are drawn in Sect. [7] 

2. Numerical Model 

The numerical model applied in our study consists of two parts. 
In a first step we simulate the hydrodynamics of the explosion 
process. Here, the description of the nuclear processes is very 
coarse. With the information gained from tracer particles ad- 
vected in this simulation we perform a nucleosynthetic post- 
processing as a second step. This enables us to infer the pro- 
duction of the individual isotopes. Both methods will be briefly 
described in the following. 

2.1. Explosion dynamics 

The deflagration model of thermonuclear supernova explosions 
as outlined in SectHlwas implem ented in a numerical scheme 
by Reine cke et alJ dl999a[l2002al) . We refer to these works for 
the details of the applied techniques and will only mention the 
basic aspects here. 

The major problem of SN la simulations is the vast range 
of relevant scales. The thickness of the flame is tiny compared 
with the dimensions of the WD star and the turbulent cascade 
interacts with the flame down to the so-called Gibson scale 
where the turbulent velocity fluctuations become comparable 
with the laminar flame speed. Neither the internal flame struc- 
ture nor the Gibson scale can be resolved in multidimensional 
simulations in the foreseeable future and thus the flame propa- 
gation and turbulence effects have to be adequately modeled in 
numerical simulations. 

Seen from the size of the WD star, it is well justified to re- 
gard the unresolved flame as a discontinuity separating the fuel 
from the ashes. Then the description of flame propagation has 
to track this interface and a technique well-suited for this pur - 
pose is the so-called level set method dOsher & Sethianl [l988). 
It is widely used in simulations of combustion problems in en- 
gineering. In this technique, the flame front is associated with 
the zero level set of a scalar field G. For numerical reasons, 
G is chosen to be a signed distance function with respect to 
the flame front. To model the flame propagation we evolve the 
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G-field according to the scheme described by Rein ecke et alJ 
d!999bh . 

In this scheme the effective flame velocity has to be pro- 
vided. To this end, the notion is essential that t urbulen t combus- 
tion proceeds in different regimes (e.g. Peters, 2000). For most 
parts of the supernova explosion the so-called flamelet regime 
applies, where the flame as a whole is wrinkled by turbulence. 
Here, the flame propagation is known to decouple from the 
microphysics of the burning proce ss and to be determined by 
the turbulent motions exclusively (Damk 5hlerl ll94ol) . These, 
however, are derived from a s ubgrid scale model implemente d 
first in SN la simulations bv iNiemever & Hillebrandll {1995). 
It describes the effects turbulence on unresolved scales. In this 
sense our model can be regarded as a Large Eddy Simulation 
(LES) well-known from computational fluid dynamics. Since 
flame propagation is modeled in our simulations, supplemen- 
tary simulations of the physical processes on small scales have 
to be provided that ens ure the validity of t he underlying as- 
sumptions. In this spirit Ropk e et al.1 J20031 l2004albh showed 
that flame propagation proceeds in a stabilized way below the 
Gibson scale. 

The hydrodynamics is modeled based on the 
PROMETHEUS implementati on JFrvxell et all 1 19891) o f 
the piecewise parabolic method (Colella & Woodward, 1984). 
The equation of state of the WD material comprises contri- 
butions from a variably degenerate and relativistic electron 
gas, ions following the Boltzmann statistics, a photon gas and 
eventually electron/positron pairs. 

The correct way to incorporate the nuclear burning would 
require a full reaction network. However, due to the restricted 
computational resources only a very simplified description 
of the nucleosynthesis is possible concurrent with the explo- 
sion simulatio n. Our imple mentation follows the approach sug- 
gested by IReinecke et al who include five species, 
viz. a-particles, 12 C, 16 0, "Mg" as a representative of interme- 
diate mass elements and "Ni" representing iron group nuclei 1 . 
The fuel is assumed to be a mixture of carbon and oxygen. 
At the initial high densities burning proceeds to nuclear sta- 
tistical equilibrium (NSE) composed of a-particles and "Ni". 
Depending on temperature and density in the ashes, the NSE 
composition chang es, which has significa nt impact on the ex- 
plosion dynamics ( IReinecke et all l2002al) . Once the fuel den- 
sity drops below 5.25 x 10 7 g cnT 3 due to the expansion of the 
WD, burning is assumed to terminate at intermediate mass el- 
ements. Below 1 x 10 7 gcirT 3 burning is switched off, since 
it is then expected to leave the flamelet regime and to enter 
the so-called distributed burning regime. Here turbulence pen- 
etrates the internal structure of the flame. This effect is ignored 
in the p resent study but was addressed by Ronke & Hilleb randl 
(l2005hh. 

In order to achieve a more detailed analysis of the nucle- 
osynthetic yields of the simulated supernova explosion we ad- 
vect tracer particles with the fluid motions recording tempera- 

1 In the following we set "Ni" and "Mg" in quotes when we refer to 
the iron group elements and intermediate mass elements followed in 
the explosion hydro-simulations. This is done to avoid confusion with 
the results of the nuclear postprocessing. 



Table 1. Nuclear reaction network (note that the elements be- 
low arsenic are irrelevant for SNe la). 
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ture, density, and internal energy as a function of time. These 
data then serve as input for a nucleosynthetic postprocessing. 

2.2. Nuclear postprocessing 

The nuclear postprocessing determines the nucleosynthetic 
yields of the explosion models a posteriori from the data 
recorded by the tracer particles. The applied m ethod is simi- 
lar to that described bv lThielemann et all ij 1986]) (there labeled 
as method (a) simple postprocessing). Its applicatio n to SN e la 
explosions is discussed in detail by Travag lio et alJ (12004). 

The employed nuclear reaction network code was kindly 
provided by F.-K. Thielemann. It comprises 384 isotopes 
which are listed in Table ^ an d takes into account /J-decays, 
electron captures, photo-disintegrations, two-body reactions, 
and three-body r eactions. A detailed des cripti on of the net- 
work i s give n by IThielem ann et al.l i 19961) and llwam oto et al! 
1 1999I) . As iBrachwitz et alJ J2000l) ~and IThielemann et al. 
(2003) discussed previously, the new electron capture and 
5-decay rates by lLanganke & Martmez-Pinedol y000) and 
Mart mez-Pinedo et alJ (l2000) are included in the network. 

Since the description of the nuclear reactions in the hydro- 
dynamic explosion simulation is coarse and Y e is assumed to be 
constant at a value of 0.5, the internal energy recorded by the 
tracer particles is employed to calculate a realist ic temperature 
form a high-temperature equation of state jTimmes & Swestvl 
2000) combi ned with an improved nuclear reaction network 
(cf. lTravaglio et alll2004l) . 
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The nucleosynthesis is calculated separately for each tracer 
particle. To level out variations in the data from the hydro- 
dynamic simulation, the minimal temperature is set to 10 9 K. 
This measure guarantees stability of the nuclear reaction net- 
work code. Subsequently, the maximum temperature r max is 
checked. If it does not exceed 2 x 10 9 K, then the correspond- 
ing material is treated as unprocessed. This approach is justified 
since the fuel consists only of 12 C, 16 0, and 22 Ne, which below 
2 x 10 9 K will burn hydrostatically, not significantly contribut- 
ing to the nucleosynthetic yields over the simulated period of 
time. 

For tracers with T mdx > 2 x 10 9 K the following procedure 
is applied: 

1 . Nuclear statistical equilibrium (NSE) is assumed if the tem- 
perature of the current time step f, is larger than 6 x 10 9 K, 
i.e. the strong reactions can be neglected and only the 
"weak" nuclear network is applied updating Y e . Otherwise 
the full reaction network is employed. 

2. Temperature and density are interpolated for the sample 
point at f, + i. If these variables change for more than 5% 
in the interval [/,, fj+i], the time step is halved. 

3. The network is solved for f, + i . If a relative accuracy of 10~ 5 
cannot be reached in a limited number of steps, the time 
step is again halved and we resume with point |2] of the 
scheme. If this measure fails, the tracer is ignored in the 
final result. Fortunately the number of such cases could be 
drastically reduced to at most one out of [27] 3 . When reach- 
ing NSE the new abundances are calculated for the updated 
Y e at t M . 

4. If the abundance of an isotope drops below 10~ 25 , it is set 
to zero. 

3. Parameter space 

The initial parameters we explore in our study (the carbon mass 
fraction X( l2 C), the central density p c , and the metallicity Z of 
the WD at ignition) are treated as independent. This allows to 
disentangle the effects of the individial parameters on the ex- 
plosion process. Nontheless, the parameter space is chosen in 
agreement with values suggested by stellar evolution, as de- 
scribed below. 

Different values for the central density of the WD and the 
carbon-to-oxygen ratio of its material are applied in the explo- 
sion model itself. Contrary to that, we vary the metallicity only 
in the nucleosynthesis postprocessing. The nomenclature of the 
models is given in Table[2] 

3.1. Variation of the carbon mass fraction 

The origin of t he diversity in the c arbon mass fraction has 
been studied bv llJmeda et alJ Jl999bll by numerically evolving 
the corresponding binary systems with 3-9 M WD progen- 
itor stars. They found it to depend on the metallicity and the 
zero-age main sequence (ZAMS) mass of the WD progenitor, 
as well as on the mass of the companion star. These in turn 
determine the mass of the WD, Mwd.o, just prior to the onset 
of accretion. The main outcome of the survey was that X( 12 C) 



Table 2. Model parameters. 
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in the core of the WD decreases with increasing Mwd.o and 
that the direct dependence of X( l2 C) on the metallicity is small 
although the correlation between the ZAMS m ass and Mwd.o 
depends sensitively on it dUmeda et all ll999a). Taking into ac- 
count the conditio ns ensuring that the W D will accrete mass 
until reaching Mrt JUmeda et al J (1199 9b) infer that X( 12 C) may 
vary in the range from ~0.36 to ~0.5. These values apply only 
to the convective core of the WD. The accreted material is as- 
sumed to be processed to a C/O ratio of ~1, leading to a gradi- 
ent of the carbon mass fraction inside the WD. This effect will 
be ignored in our model, were we postulate a uniform C/O ratio 
throughout the entire star employing values of 0.30, 0.46, and 
0.62 for X( 12 C) (cf. TableQ. 

3.2. Variation of the central density 

The variation of the central density in SN la progenitors just be- 
fore the ignition of the flame is even more difficult to constrain. 
At least two effects determine the value of p c . 

The first i s the accretion history of the binary system (see 
Lan ger et all l2000l for a detailed study of the accretion pro- 
cess). There seems to be only a narrow window in the range of 
possible accretion rates M in which carbon can be ignited cen- 
trally, avoiding off-center ignition s and gravitational coll apse 
due to high electron-capture rates. Nomoto & Iben ( 1985| re- 
port on two centrally ignited models with p r — 1 .7 x 10 9 g cm -3 
and p c — 5.2 x 10 9 g cirT 3 , respectively, and lBravo et all ill 9961) 
find models in the range 1.8 x 10 9 gcirr 3 < p c < 6.3 x 
10 9 g cm~ 3 . However, the exact range of that window is uncer- 
tain and depends additionally on the white dwarf mass and tem- 
perature. Initially cooler WDs are shifted to rather high central 
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densities in the range 6 x 10 9 gcm 3 < p c < 1.3 x 10 I0 gcm 3 
faernanz et allll988l) . 

The second effect is the establishment of the thermal struc- 
ture of the WD. Cooling due to plasmon neutrino losses 
and neutrino b r emsstr ahlung have to be taken into account 
Jlwamoto et all Il999l) . and a (most uncer tain) contribution 
may come from the convective Urea process JPaczvhsklll972| 
Bar kat & Wheelerl Il990t iMochkovitchl Il996l iLesaffreet all 



Brav o et alJ 41,993) calculate models for central densities 
at ignition of 2. 5 x 10 9 gem' 3 4.0 x 10 9 gcirT 3 , and 8.0 x 
10 g cm " 3 , and llwamoto et all {l999) use values of 1.37 x 
10 9 gcrrT 3 and 2.12 x 10 9 gcrrT 3 . We assume central densi- 
ties of 1.0 x 10 9 gcirT 3 and 2.6 x 10 9 gcirT 3 (see Table E}. 
Unfortunately it is not yet possible to apply higher central den- 
sities, since then electron captures will become dynamically 
important. Although electron captures are correctly treated in 
the nuclear postprocessing, this effect is not implemented in 
the current explosion models. 

3.3. Variation of the metallicity 

Our ignorance concerning realistic progenitor evolution is 
evident in the approach to prescribe the metallicity Z 
of the progenitor independent of the other parameters- 
Detailed stellar evolution c alculations (e.g. lUmeda et alJ 
ll999bHDomfnguez et alJl200lh have shown that it strongly in- 
fluences the progenitor's central density and also the C/O ratio. 
Nevertheless, in the spirit of our exploration of possible effects 
in the explosion models, we set aside a realistic progenitor de- 
scription and treat the metallicity as an independent parameter. 
A direct effect of the metallicity of the WD's progenitor is the 
14 N abundance after the CNO burning phase. In Helium burn- 
ing it is then converted mostly to 22 Ne. For simplicity, we as- 
sume a uniform distribution of this 22 Ne, which is only justified 
in regions mixed by pre-ignition convection. 

An analytic estimation of the effect of the metallicity on the 
56 Ni production was given bv lTimmes et alJ J2003h . They sug- 
gest a variation of Z ranging from 1 /3 to 3 times solar, based 
on observations of field white dwarfs. Following this sugges- 
tion we vary the 22 Ne abundance in our models to simulate a 
variation in metallicity Z. In particular we explore Z (corre- 
sponding to X( 22 Ne) = 0.025), 0.5Z o , and 3Z G (cf. Table©. 

4. Simulation setup 

A rather large number of simulations is required to cover the 
parameter space we aim at. We therefore have to minimize the 
computational expenses by applying a simple setup for the in- 
dividual models. Our calculations span only one spatial octant 
and assu me mirror symmetry to the oth er octants. Full-star sim- 
ulations ( Ropke & Hillebrandt, 2005a) have shown that this ap- 
proach does not miss large-scale flame features and thus - al- 
though being a simplification - does not restrict the validity of 
the model. The simulations were set up on a Cartesian com- 
putational grid that was equally spaced in the inner regions. 
To capture the expansion of the WD, the outer g rid cells were 
widened exponentially. Recently, Ropke (2005) showed that 



with a comoving computational grid the evolution can be fol- 
lowed to homologous expansion. This, however, is not applied 
in the present models. 

The resolution of the individual runs was rather low - the 
computational domain was divided in [256] 3 grid cells corre- 
sponding to a central grid resolution of 10 6 cm. In each di- 
rection the grid length in the outer 35 zones was increased 
subsequently by a fac tor of 1.15. As was pointed out by 
IReinecke et alJ J2002cl) the chosen resolution still guarantees 
the explosion characteristics to be numerically converged (pos- 
sibly with the exception of the latest stages of the burning 
where intermediate mass elements are produced). However, 
with this resolution it is not possible to set up reasonable multi- 
point ignition scenarios, since only a very small number of 
seed-bubbles could be r esolved . This is certainly a drawback 
because lkeinecke et alJ J2002bl) showed that such models give 
rise to more vigorous explosions. We restrict our simulati ons to 
the cent rally ignited model c3_3d_256 model o flReinecke et alJ 
(2002c) in which the spherical initial flame geometry is per- 
turbed with three toroidal rings (see the upper left panel of 
Fig-Q- Note that we initially incinerate the same volume in all 
models, which does not correspond to the same mass for dif- 
ferent central densities. This ensures the same initial numerical 
resolution of the flame front. 

For the construction of the WD near the Chandrasekhar 
mass we follow the procedure described by Reinecj3 d200ll) . 
We assume a cold isothermal WD of a temperature Tq = 
5 x 10 5 K. With the chosen values for the carbon mass frac- 
tion of the material and the central density we integrate the 
equations of hydrostatic equilibrium using the equation of state 
described in Sect. 12.11 Depending on the central densities and 
compositions the masses of the resulting WDs vary slightly: for 
p c = 1.0X 10 9 gcirT 3 andp c = 2.6 x 10 9 gcirT 3 the WD masses 
amount to 1.367 M and 1 .403 M , respectively. As tested by 
IReinecke! feoOl). the construction procedure guarantees stabil- 
ity of the WD over a time longer than simulated. 

The [«trace] 3 tracer particles are distributed in an n traC e X 
"trace x "trace equidistant grid in the integrated mass Mo(r), the 
azimuthal angle 0, and cos 6, so that each particle represents 
the same amount of mass. In order to improve the tracer par- 
ticle statistics, a random offset to the coordinates was applied. 
This offset was chosen small enough to keep the tracer parti- 
cles in their individual mass cells. The values of the density, the 
temperature and the internal energy at the tracer particle's loca- 
tion and its coordinates were recorded every ~1 ms. This allows 
for an accurate reconstruction of the trajectories as well as the 
final velocities and the thermodynamical data. In the models 
presented in the following we set « trace = 27. To test the rep- 
resentation of the model in the tracer particles in cases of low 
central densities, this number was increased to 35 in test calcu- 
lations, as will be discussed below. 

5. Explosion models 

The explosion simulation for the exemplary case of model 
2_2_X (the metallicity does not affect the explosion dynam- 
ics in our implementation) at four different times is illustrated 
in Fig. ^ The isosurface indicating the position of the flame 
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Fig. 1. Time evolution of the burning front for model 2_2_X. 



front is rendered from the zero level set of the scalar field G. 
The computational grid plotted in these snapshots visualizes 
our setup with uniform grid cells in the inner region and an 
exponential growth of the grid spacing further out. Our ini- 
tial flame configuration is shown in the upper left snapshot of 
Fig. \l\ In the subsequent snapshots the growth of instabilities 
and an increasing wrinkling of the flame front are visible. Once 
the flame enters the exponentially growing part of the grid, the 
resolution of flame features becomes coarser. However, at this 
stage the expansion of the WD decreases the density of the fuel 
to values where burning has largely ceased in our model. Thus 
the coarse flame resolution in late stages of the simulation does 
not affect the results. 




Fig. |2 shows the total energy production of our models. 
Due to the simple setup, all explosions are weak, but trends 
can clearly be identified. The energy releases of the different 
models are listed in Table [3] which also provides the masses 
of produced iron group elements ("Ni") and intermediate mass 
elements ("Mg"). In Figs. [3] and |4] the energy generation rates 
and the evolution of the turbulent energies in our models are 
plotted, respectively. 



-6 I i i i i I i i i i I i i i i I i i i i 

0.0 0.5 1.0 1.5 2.0 

t[s] 

Fig. 2. Total energies in models (a) 2_3_X, (b) 2_2_X, (c) 
2_1_X, (d) 1_3_X, (e) 1_2_X, and (f) 1_1_X. 
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Table 3. Results of explosion models: produced masses of iron 
group elements ("Ni") and intermediate mass elements ("Mg"), 
nuclear energy release, and total energy at the end of the simu- 
lations. 



model 


M("Ni") 


M("Mg") 




Em 




[M ] 


[M Q ] 


[10 50 erg] 


[10 50 erg] 


lJJt 


0.3944 


0.2067 


6.974 


2.714 


1_2_X 


0.3867 


0.2081 


7.445 


3.140 


1_3_X 


0.3757 


0.2144 


7.870 


3.563 


2_1_X 


0.5178 


0.1874 


8.851 


3.772 


2_2_X 


0.5165 


0.1859 


9.461 


4.412 


2_3_X 


0.5104 


0.1822 


9.966 


4.909 



20 




0.0 0.5 1.0 1.5 2.0 

f[s] 

Fig. 3. Energy generation rates in models (a) 2_3_X, (b) 2_2_X, 
(c) 2_1_X, (d) 1_3_X, (e) 1_2_X, and (f) 1_1_X. 



5.1. Variation of the progenitor's C/O ratio 

The effects of a variation of the progenitor's carbon-to-oxygen 
ratio on the SN la explosio n models have been described by 
iRopke & Hillebranddd2004T) . We extend the discussion here. 

Considering the explosion energetics first, Fig.|5]shows that 
a higher carbon mass fraction leads to an increased energy pro- 
duction for fixed central densities. Values are given in Table|3] 
For both central densities the nuclear energy releases of the 
models increase by 12% (~27% in the total energies) changing 
X( l2 C) from 0.30 to 0.62. The observed trend is not surpris- 
ing and can easily be explained by the burning process. The 
predominant effect is certainly the difference in the mean bind- 
ing energy of the fuel. A higher carbon mass fraction increases 
the total energy generation for the simple reason that the bind- 
ing energy of 12 C is lower than that of lfi O so that it releases 
more energy by fusion to iron group elements. A minor ef- 
fect coul d be that the laminar bur ning velocity increases with 
X( l2 C) dTimmes & Woo slev. 1992). This, however, is negligi- 
ble in our models, since already after a few time steps the flame 
propagation is completely determined by the turbulent flame 
speed. 



4r 




0.0 0.5 1.0 1.5 2.0 

f[s] 

Fig. 4. Turbulent energies in models (a) 2_3_X, (b) 2_2_X, (c) 
2_1_X, (d) 1_3_X, (e) 1_2_X, and (f) 1_1_X. 

It is noteworthy that the evolution of the energetics in the 
model does not show a strong temporal shift. The energy gen- 
eration rate peaks at comparable times for the models with dif- 
ferent carbon mass fractions (cf. Fig. [3] for the temporal evo- 
lution of the energy generation rate 2 ). Fig. |2 reveals that the 
total energies of our models are very similar for the largest part 
of energy generation and differ o nly in the late pha ses. In this 
point our findings disagree with Khokhlov (2000). Although 
he speculates that a decreasing X( l2 C) would result in weaker 
explosions, he claims that the reason is a delay of the develop- 
ment of the buoyancy instabilities, which seems to be only a 
minor effect in our simulations. The reason for the difference 
in the interp retation of the results is possibly that the models of 
iKhokhlovllEoOol) were apparently not evolved beyond the max- 
imum of energy generation so that it is difficult to distinguish 
between a delay and an overall lower energy production. 

Contrary to the explosion energetics, the produced masses 
of iron group elements are u nexpected. The working hypoth- 
esis by Umed a et alJ (11999a) predicting a larger production of 
56 Ni for larger carbon mass fractions was based on the specula- 
tion that the resulting increased energy release would enhance 
buoyancy effects and thus accelerate the turbulent flame prop- 
agation. Consequently, more material would be burnt at high 
velocities producin g larger fract i ons of iron group elements. 
As emphasized by Um eda et alJ dl999al) . this hypothesis can 
only be tested in multi-dimensional simulations which treat the 
turbulent flame velocity in an unparametrized way. This is pro- 
vided by our approach, but surprisingly our models do not sup- 
port the hypothesis. The energies in our models evolve simi- 
lar in the first part of the explosion and no enhanced turbulent 
flame propagation is visible regardless of the carbon mass frac- 
tions. The similar temporal evolutions of the energetics in our 
models correspond to a striking similarity in the flame mor- 
phology. Fig.[5]illustrates the flame front in models with differ- 

2 Note, that the peak at t = s is caused by our setup in which the 
initial flame is initialized by instantly incinerating the material behind 
it. 
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model 2_2_X 




Fig. 5. Flame surface of models with different carbon mass fraction at t — 1.0 s. 



1.2 



model 2_1_X 
model 2_2_X 
model 2 3 X 




Fig. 6. Mean effective gravitational accelerations experienced 
by the flame fronts for models with different C/O ratios. 



ent C/O ratios at t — 1.0 s. The extent of the burnt volumes is 
comparable. The similarities in the first part of the energy gen- 
eration as well as in the flame morphologies and flame propaga- 
tions indicate that the large scale buoyancy effects, which feed 
energy into the turbulent cascade and thereby drive the flame 
propagation, are comparable in models with different C/O ra- 
tio. The buoyant velocities can be estimated from the relation 



^buoy 



y/AtgL 



(1) 



for a single non-burning rising bub ble of size L subj ect to 
a gravitational acceleration g (Davies & Tavlor, 1950). The 
Atwood number A t characterizes the contrast between the den- 



sity inside the bubble (pi) and outside it (p Q ): 

bo-Pil 



At 



Po +Pi 



(2) 



In a supernova explosion, the situation is, of course, much more 
complex since bubbles burn and will merge. Nevertheless, it is 
clear that the effective gravitational acceleration (Atg) deter- 
mines the flame propagation velocity - not only directly on the 
largest scales but also over the mechanism of the interaction 
of the flame with generated turbulence. This effect can only be 
revealed in multidimensional calculations, as presented here. 
Fig-EDshows the mean effective gravitational acceleration (At g) 
experienced by the flame front. Only minor differences are vis- 
ible here. The data point at t — 0.0 s is unphysical, since the 
material behind the flame had not been burnt at this instant and 
thus there is no density contrast over the flame yet. With tem- 
poral evolution there is a competition between a rapidly de- 
creasing gravitational acceleration due to the expansion of the 
star in the explosion and an increasing de nsity contrast over the 
flame along with lower fuel densities (cf. Tim mes & Wooslevl 
119921) . As seen from Fig. [6] finally the decreasing gravitational 
acceleration dominates this competition. 

Because of the very similar evolutions of the large-scale 
buoyancy effects, there are little differences in the evolutions 
of the turbulent energies in models with varying carbon mass 
fractions (see Fig.|4f). 

If the larger energy from burning carbon-rich fuel was di- 
rectly converted to work expanding the ashes, buoyancy effects 
and an acceleration of the turbulent flame propagation should 
be observable in our simulations. The only way to bypass these 
effects is that the energy is buffered in a larger fraction of a- 
particles present in the NSE. This fraction indeed increases 
with higher temperatures and the consequences are twofold. 

3 The values of £ turb are not significant at late times since those 
are derived from the subgrid energy which depends on the length of 
the grid cells. In the outer regions which the flame enters at late times, 
those are elongated and therefore E tmh rises again after reaching a peak 
at t ~ 0.65 s and t ~ 1.05 s for p c = 2.6 x 10 9 gem -3 and p c = 1.0 X 
10 9 gcirT 3 , respectively. For a uniform grid it would be expected to 
monotonically decrease after these peaks (cf. Rooke, 2005). 
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Fig. 7. Fuel densities at the mean flame locations as a function 
of the mass of "ashes" behind the flame front for models with 
different carbon mass fraction of the progenitor material. 




Fig. 9. Mean effective gravitational accelerations experienced 
by the flame fronts for models with different central densities. 
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Fig. 8. Temporal evolution of the chemical composition in 
models with different carbon mass fraction of the progenitor 
material. 



First, the binding energy of the ashes is lowered and less en- 
ergy is released from thermonuclear burning. Second, distribut- 
ing the energy on an increased number of particles in the ashes 
decreases their temperature. Both effects lead to an increase in 
the density of the ashes which suppresses the buoyancy effects. 
Hence the turbulent flame propagation velocity in carbon-rich 
fuel models is lowered to values comparable to those found in 
oxygen-rich fuel simulations. As a consequence similar masses 
of fuel are burnt at particular fuel densities. To corroborate this, 



we plot the fuel density at the average flame front location over 
the burnt mass in Fig.0 

Fig- 13 proves that the effect of energy buffering in the a- 
particles indeed applies for our models, here shown for the 
models with a central density of 2.6 g cirT 3 . Between 0.2 s and 
0.9 s the ashes contain significant amounts of c-particles. The 
maximum mass fraction of c-particles increases by about 20% 
when changing the carbon mass fraction from 0.30 to 0.62. 
This effect is capable of compensating the differences in the 
fuel binding energies according to the following estimate. At 
t ~ 0.6 s (the energy generation rate peaks here, cf. Fig. O the 
ashes in the model 2_2_X contain about 21% or-particles and 
79% "Ni". If there was no change in the amount of a-particles 
along with changing the C/O ration in the models, the nuclear 
energy release would increase for about 22% changing X( l2 C) 
from 0.30 to 0.62. In contrast, taking into account the observed 
20% increase in the a-particles in the ashes, the nuclear energy 
release difference reduces to ~5%. 

This self-regulation mechanism has an important conse- 
quence. Since it suppresses increased buoyancy effects which 
would otherwise arise with larger carbon mass fractions in the 
fuel, similar amounts of fuel are consumed by the flame at 
stages where burning terminates in NSE. Therefore the amount 
of produced iron group elements is nearly kept constant for dif- 
ferent C/O ratios in the fuel. 

The a-particles buffer the energy only temporarily. With 
further expansion and cooling of the WD they are converted to 
iron group elements ("Ni") releasing the stored energy. This is 
the reason why the energies in the models diverge at later times. 
Then, however, the fuel density has dropped to values where 
burning terminates in intermediate mass elements and hence 
the synthesis of iron group elements is unaffected. Therefore 
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the models with different carbon-to-oxygen ratios produce sim- 
ilar amounts of iron group elements. Interestingly, we find even 
a slight decrease in the production of iron group elements for 
an increasing carbon mass fraction of the model. The same 
holds for the intermediate mass nuclei in the high central den- 
sity models while the trend is reversed in the models with lower 
central densities (cf. Table[3j. 

5.2. Variation of the central density 

For higher central densities at ignition, the explosion turns out 
to be more vigorous for a fixed carbon mass fraction in the WD 
material (cf. Fig. 13 Table[3}. Here, the nuclear energy releases 
differ about 26% and the total energies vary about 34%. 

As for an increased carbon mass fraction, a higher 
density of the fuel acc elerates the laminar burning 
jTimmes & Wooslevl [f992). Again, this has little impact 
on our models since burning proceeds laminar only in the first 
few time steps. Two other effects are more significant here. 

First, for the higher central density, obviously more mate- 
rial is present at sufficiently high densities so that it can poten- 
tially be burnt to iron group elements. 

Second, for higher central densities the effective gravita- 
tional acceleration (At g) experienced by the flame is higher in 
the first ~0.9 s (cf. Fig. |9jl- This increases the development of 
flame structures resulting from the non-linear stage of the buoy- 
ancy instability. As a result, the turbulent cascade will develop 
more quickly (cf. Fig. |4j and the turbulence-induced boost of 
of the effective flame propagation velocity sets in earlier. This 
is shown in Fig.^3 where snapshots of the flame evolutions at 
fixed times for models with different central densities at igni- 
tion are given. 

Consequently, the production of iron group elements in- 
creases with higher central densities, while the amount of in- 
termediate mass elements decreases. 



6. Nucleosynthesis 

After the nucleosynthesis postprocessing the abundances of the 
individual isotopes in the ashes are unveiled. Additionally, the 
parameter of the progenitor's metallicity comes into play here 
by assuming a certain fraction of the material to be composed 
of 22 Ne. 

The lightcurves of SNe la are powered by the radioac- 
tive decay of 56 Ni and 56 Co. Therefore their peak luminosi- 
ties are determined by the nucleosynthetic yields of the ex- 
plosions rather than the energetics. Consequently, the nucle- 
osynthetic postprocessing of explosion models can shed some 
light on the observed diversity of SN la events. 56 Co decay is 
slow and thus the peak luminosity is a function of the produced 
amount of 56 Ni. A compilation of the 56 Ni masses derived from 
all our models by nucleosynthetic postprocessing can be found 
in TableH 

Although we will focus here on the effect of initial param- 
eters on the production of 56 Ni, we will first discuss the overall 
picture of the nucleosynthesis yields. 



6.1. The final yields 

The freeze-out masses after completion of the /^-decays are 
plotted in Figs. ^2 to Elf or different models. Here the usual 
normalization to the solar abundances and the 56 Fe mass frac- 
tion was applied. Values are given in Table[6]in the online ma- 
terial. 

Fig. II llshows a comparison between models with different 
carbon mass fractions for fixed p c = 2.6 x 10 9 gcrrT 3 and for 
solar metallicity. Obviously, the carbon mass fraction has only 
little effect on the final abundances. Though some variation is 
visible for the intermediate mass elements (Mg to Ca), there is 
practically no difference in the iron group yields for the differ- 
ent models. This is expected from the analysis of the explosion 
process in the previous section. Due to the energy buffering 
in the a-particles, burning to NSE consumes almost identical 
masses of fuel, while the recombination of the a-particles at 
the end of complete NSE burning leads to an additional energy 
release that varies with the C/O ratio. Therefore the incomplete 
burning in the models that follows burning to NSE proceeds 
differently in the various models. 

The variations in the final yields due to different central 
densities are illustrated in Fig.^] Here, the models X_2_2 are 
plotted, i.e. the C/O ratio is fixed to 0.81 and the metallicity is 
solar. The model with lower central density produces more in- 
termediate mass elements, but the variations are small. In con- 
trast, for higher central densities, there is a visible increase in 
the abundances of iron group isotopes, viz. titanium, vanadium, 
chromium, manganese, iron, cobalt, and nickel. The two effects 
that contribute to an increased mass consumption in complete 
NSE burning with higher p c were discussed in Sect. 15.21 The 
resulting final yields are a natural consequence of these effects. 

Changes in the progenitor's metallicity resulting in differ- 
ent abundances of 22 Ne in the WD material have a large im- 
pact on the final yields. To illustrate this influence, we consider 
the models 2_2_X for 0.5 Z , 1.0 Zq, and 3.0 Z (cf. Fig. P. 
The variation of the 22 Ne abundance is obvious and caused by 
the representation of the progenitor's metallicity in the different 
mass fractions of that isotope in our simulations. The produc- 
tion of chromium, manganese and iron isotopes is increased for 
higher metallicity, especially for isotopes with two more neu- 
trons than the symmetric nuclei ( 54 Fe, 58 Ni). An exception is 
56 Fe which was used to normalize the abundances. This trend 
holds analogously for intermediate mass elements. In particu- 
lar, one observes a higher abundance of 26 Mg, 30 Si, 34 S, 38 Ar 
and 42 Ca with increased metallicity. Comparing the models 
2_2_2 and 2_2_3 the change is by a factor of 1 1 for 26 Mg and a 
factor of approximately 3 for the other isotopes (cf. Table[6]in 
the online material). The other models not present in the table 
give similar factors for identical metallicities. The increase of 
neutron-rich isotopes is caused by the fact that a higher pro- 
genitor's metallicity results in an increased 22 Ne mass fraction, 
which serves as a source of a neutron excess. 
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Fig. 10. Flame surface of models with different central densities prior to ignition at t — 1.0 s. 



model 


p c [10 9 gcnr 3 ] 


X( [b O) 


X( [2 C) 


X( 22 Ne) 


Z[Z ] 


M( 56 Ni)[M ] 


UJ 


1.0 


0.7 


0.2875 


0.0125 


0.5 


0.2982 


1_1_2 


1.0 


0.7 


0.275 


0.025 


1.0 


0.2876 


UJ 


1.0 


0.7 


0.225 


0.075 


3.0 


0.2450 


UJ 


1.0 


0.54 


0.4475 


0.0125 


0.5 


0.2966 


1_2_2 


1.0 


0.54 


0.435 


0.025 


1.0 


0.2860 


1J2_3 


1.0 


0.54 


0.385 


0.075 


3.0 


0.2444 


UJ 


1.0 


0.38 


0.6075 


0.0125 


0.5 


0.2907 


1_3_2 


1.0 


0.38 


0.595 


0.025 


1.0 


0.2805 


1_3_3 


1.0 


0.38 


0.545 


0.075 


3.0 


0.2403 


2_1_1 


2.6 


0.7 


0.2875 


0.0125 


0.5 


0.3115 


2_1_2 


2.6 


0.7 


0.275 


0.025 


1.0 


0.2999 


2_1_3 


2.6 


0.7 


0.225 


0.075 


3.0 


0.2544 


2_2_1 


2.6 


0.54 


0.4475 


0.0125 


0.5 


0.3163 


2_2_2 


2.6 


0.54 


0.435 


0.025 


1.0 


0.3046 


2_2_3 


2.6 


0.54 


0.385 


0.075 


3.0 


0.2592 


2_3_1 


2.6 


0.38 


0.6075 


0.0125 


0.5 


0.3174 


2_3_2 


2.6 


0.38 


0.595 


0.025 


1.0 


0.3065 


2_3_3 


2.6 


0.38 


0.545 


0.075 


3.0 


0.2608 



Table 4. 56 Ni masses synthesized according to the nucleosyntehsis postprocessing. 



6.2. Impact of the CIO ratio on the 56 Ni mass 

Contrary to the previous section, we analyze here the nucle- 
osynthesis yields right after the explosion. The production of 
radioactive species is given in Table[7]in the online material. 

Fig. El shows the 56 Ni production for the models in de- 
pendence of the carbon mass fraction of the progenitor. The 
central densities are fixed top f = 1.0 x 10 9 g cm and p r = 
2.6 x 10 9 gcrrT 3 , respectively. The metallicities of the models 
shown here are set to Z = Z . We note only minor changes in 
the 56 Ni masses (about 2%) for both central densities, which is 
not surprising given the small variations in the flame morphol- 
ogy and advancement discussed in Sect. 15. 21 

It should be noted in Fig. [21 that the trend of 56 Ni pro- 
duction has opposite directions for different central densities. 



While this feature is in accordance with the total production 
of iron group elements in the explosion models for low cen- 
tral densities, it is reversed for the high central density case (cf. 
Tabled- In order to check whether an under-representation of 
NSE-material in tracer particles in the low density case was the 
origin, we recalculated these models with the number of tracers 
increased to 35 3 . The trend of decreasing 56 Ni production with 
higher carbon mass fraction was weaker, but had still the same 
direction. Since the variations are at the percent level, it is be- 
yond the accuracy of our models to judge whether the trend is 
of physical nature or an artifact of our simulation. 

The result of the 56 Ni production in the explosion phase 
being largely independ ent of the carbon mass fraction supports 
the conjecture of Ropke & Hillebrandtl d2004l) that the peak lu- 
minosity of SNe la will be only marginally affected by the 
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Fig. 11. Final abundances for models with different C/O ratios. 
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Fig. 14. 56 Ni production depending on the carbon mass fraction 
for models with different central densities and solar metallicity. 



carbon-to-oxygen ratio of the progenition WD star. This con- 
jecture was only based on the cumulative production of all iron 
group elements and is now confirmed by the derivation of the 
exact amounts of 56 Ni via nucleosynthetic postprocessing. 



6.3. Impact of the central density on the 56 Ni mass 

For fixed C/O ratios of 0.81 and solar metallicities our mod- 
els produce 0.286 M Q of 56 Ni for p c = 1.0 x 10 9 gem and 
0.305 M Q of 56 Ni forp c = 2.6 x 10 9 g cirT 3 , i.e. from the lower 
to the higher central density the 56 Ni production increases for 
7%. These changes go along with the higher overall production 
of iron group nuclei at higher central densities (cf. Table |3}. 
The reasons for this effect have been discussed in Sect. 15.21 

Although somewhat larger than the changes found in the 
case of varying carbon mass fraction, the effect is still rather 
small. However, our study covers only parts of the effects re- 
sulting from changing the central densities of the models. With 
a further increasing central density, electron captures will be- 
come important and the 56 Ni production is expected to decrease 
while the total mass of iron group elements should still in- 
crease. Unfortunately, in the current study this effect could not 
be consistently modeled, but it will be addressed in forthcom- 
ing work. 

6.4. Impact of the metallicity on the 56 Ni mass 

iTimmes et alJ J2003h proposed an analytic model for the effect 
of the progenitor's metallicity on the 56 Ni production in SN la 
explosions. Their reasoning is based on the assumptions that 
most of the 56 Ni is produced between the 0.2M Q and the O.8M 
mass shells in NSE and that the Y e is constant during burning 
in that region. This is motivated by one-dimensional models. 
Furthermore, they take into account only the species with the 




Fig. 12. Final abundances for models with different central densities. 
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with Z = Z/Z . This equation is obtained from combining the 
equations 

M56 N ; 



M 3 6 Ni (Z = 0) 



: 58K - 28, 



resulting from conservation of mass and charge, 



X( 22 Ne) = 22 



X( n C) X( H N) X( 16 0) 



12 



14 



16 



(4) 



(5) 



approximating the 22 Ne abundance resulting from the metallic - 
ity of the ZAMS progenitor, and 
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Fig. 15. 56 Ni production depending on the progenitor's metal- 
licity for models with different central densities C/O ratios. 



highest mass fraction; in a first step 56 Ni and 58 Ni. Under these 
assumptions they derive a linear correlation between Z and the 
produced M( 56 Ni): 



10 „ 26 , fi 
Y e = 2^X( 22 Ne) + ^X( 56 Fe) 

+ -[l - X( 22 Ne) - X( 56 Fe)] . 



(6) 



giving the initial Y e of the white dwarf right before the explo- 
sion under the assumption of a uniform distribution of 22 Ne and 
56 Fe. When the presence of 54 Fe is take n into accoun t the factor 
0.057 in Eq. Q changes to 0.054 (cf. lTimmes etalll2003l) . 

For a comparison of this analytic prediction with our mod- 
els we set X( 56 Fe) — in Eq. (|5Jl since here the initial Y e is 
determined by 22 Ne only. We set Z( 22 Ne) = 0.025Z, where 
Z = Z/Zq. If Eq. (|5Jl is now substituted into Eq. 0, the follow- 
ing equation is obtained 



M* Ni (Z) 
M 56Ni (Z = 0) 



1 - 0.057Z, 



(3) 



gggg = l-g» 0.Q252. 

Ms6 Ni (Z = 0) 22 



(7) 
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Fig. 13. Final abundances for models with different metallicities. 



To compare this linear dependence of the 56 Ni mass on Z 
with our simulations, a linear regression following 

M 5 6 Ni (Z) = M + mZ (8) 

was applied to our data. Here Mq denotes the extrapolated value 
M( 56 Ni) at Z = 0. Values for m and M for the different models 
are given in Table[5] The coefficient of correlation is 0.9999 in 
all cases, which suggests a good agreement of our data with a 
linear dependence, but, of course, more data points would be 
desirable for a definite statement. 

To compare the slope of 0, i.e. -0.0659, with the fits to our 
data, we give values for -m/Mo in Table [5] The agreement is 
reasonable keeping in mind that Eqs. and @ were derived 
by assuming 56 Ni and 58 Ni to be the two most abundant iso- 
topes in NSE with a constant Y e . However, a significant amount 
of 56 Ni seems to be produced in regions where the assumption 
of constant Y e breaks down. 

Th us, the analytical model introduced in iTimmes et alJ 
provides an excellent explanation for the effect of the 
metallicity. Based on models diffe rent from the ones app lied 
here this was recently confirmed by Trava glk) et alJ J2005 ). 

7. Conclusions 

In the present paper the impact of several progenitor parameters 
on three-dimensional SN la explosion models has been studied 
for the first time in a systematic way. Here, we investigated 
the effects of the progenitor's central density, its carbon mass 
fraction, and its metallicity. Of course, there may be several 



other parameters that possibly affect the light curve of SNe la 
(rotation of the progenitor, morphology of the ignition spot(s), 
a delayed detonation at varying densities, asphericities etc.), 
which were not addressed in the present survey. 

A first important point to note is that our numerical imple- 
mentation as well as the underlying astrophysical model are 
evidently robust against variations of the initial conditions to a 
reasonable degree. On the one hand, the variations in the result- 
ing features are relatively small. A deviation in orders of mag- 
nitude should have been reason for concern, but all our models 
seem to be well-behaved. On the other hand, the model is not 
too robust in the sense that variations of the initial parameters 
do show effects on the results, i.e. an intrinsic variability is pre- 
served. The degrees of freedom expected for a SN la explosion 
are at least not entirely artificially suppressed in our model. 
Hen ce our mo del fulfills the req uirements 2, 3, and 4 stated by 
iHillebrandt & Niemeveil l l2000l) . 

Another point is the absolute scale of the results. Given the 
limited resources of computational time and storage space, we 
had to restrict the models to a resolution of [256] 3 grid cells 

per octant. Although suc h models reach n umerical convergence 

I 1 1 I 

in global characteristics (Reinecke et al., 2002c), it is not pos- 
sible to apply multi-spot ignition scenarios at this resolution 
which would produce more vigorous explosions. As a conse- 
quence, the explosion energy of all our models is rather low 
and the 56 Ni productio n falls short the nickel mass of a pro- 
totype SN la JContardo et alJ feOOOl) find O.41M of 56 Ni for 
SN1994D). These restrictions exclude the possibility of find- 
ing the absolute scale of effects and hence requirement 1 of 
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model 


1 1 X 


1 2 X 


1 3 X 


2 1 X 


2 2 X 


2 3 X 


— m 


0.0213 


0.0209 


0.0201 


0.0228 


0.0228 


0.0227 


M [M ] 


0.3089 


0.3070 


0.3007 


0.3228 


0.3275 


0.3290 


-m/M [M- 1 ] 


0.0690 


0.0681 


0.0668 


0.0706 


0.0696 


0.0690 



iHillebrandt & Niemeverl d2000l) is not met in the current study. 
However, there is a fair chance that models with more elabo- 
rate initial flame representations will agree better with t he ab- 
solute values of observed quantities (see e.g. iTravaelio et all 
2004). Nevertheless the present parameter study should reveal 
the correct trends of the variation of SN la properties. 

A major uncertainty lies in the range of variation of the pro- 
genitor parameters. Although we applied values that are com- 
mon in literature, our parameter space is not derived from a 
realistic stellar evolution of the progenitor. 

Keeping this in mind, the maximum variation in 56 Ni of 
about 27% found in our parameter study can be regarded as a 
strong hint that the variations of the progenitor properties taken 
into account here provide a significant contribution to the scat- 
ter in SN la luminosities. However, it seems unlikely that these 
are sufficient to explain the full range of diversities in "Branch 
normal" SNe la. Of course, more elaborate models are required 
to assess this. 

Regarding the diversity of 56 Ni production in our models 
resulting from the variation of the initial parameters, the fol- 
lowing trends were found: 

- The progenitor's carbon-to-oxygen ratio has only little im- 
pact on the amount of produced 56 Ni. This is in strong con- 
trast to the common assumption that the C/O ratio be a ma- 
jor source of luminosity variation in SN la e xplosions. The 
"working hypothesis" o fllJmeda et alJdl999al) could not be 
confirmed by our models. The reason for this effect could 
be unveiled by our three-dimensional simulations. Since 
flame propagation in the deflagration stage is mainly de- 
termined by the turbulent motions of the material, the ex- 
plosion dynamics is not altered as long as the buoyancy 
effects that generate the turbulence are comparable. This is 
given in our models at stages of iron group nuclei synthesis. 
Different energy releases resulting from differences in the 
fuel binding energies are compensated by a varying amount 
of a-particles present in the ashes. These buffer the temper- 
ature of the ashes and thus the densities are not altered sub- 
stantially ensuring the same buoyancies. Consequently, the 
explosion dynamics is similar in the stages of iron group 
element synthesis for models with different C/O ratio in the 
fuel resulting in a small variation of the produced 56 Ni of 
only about 2%. 

- The central density affects the 56 Ni production. The vari- 
ation found in our models amounts to about 7%. This is 
explained by the fact that for higher central densities more 
material is burned under conditions where iron group el- 
ements are produced. Moreover, a higher central density 
increases the mean gravitational acceleration experienced 
by the flame front and thus enhances the generation of tur- 



bulence thereby accelerating the flame propagation. Due to 
this effect even more material is processed at higher densi- 
ties where the reactions terminate in iron group elements. 
- A greater effect (assuming that our parameter space is rea- 
sonable) was found for a variation of the metallicity in the 
nuclear postprocessing. By varying the 22 Ne mass fraction 
from 0.5 to 3 times solar, a variation of the produced 56 Ni 
mass of about 20% was found. Our models were consistent 
with the analytical prediction of a linear relation b etween 
the metallicity and X( 56 Ni) bv lTimmes et alJfe003l) . 

The effects of varying C/O ratios and central densities of the 
progenitor on the supernova explosion are based on effects of 
the turbulent flame propagation and can thus only be revealed 
by three-dimensional models. 

However, we have to emphasize an important limitation of 
the results. Our analysis addresses only changes in the explo- 
sion process itself. For comparability of the simulations we as- 
sumed identical initial flame configurations. The ignition pro- 
cess, however, may be influenced b y the c arbon-to-oxygen ra- 
tio of the progenitor dWooslev et alll2004 . Since different ini- 
tial flam es can have large impact on the explosion dynamics 
(e.g.lRei necke et aHl2002l l lGam^z o et allE()03UCalder et all 
2004i lRopke & HillebrandtU2005al) . theC/O ratio may still be 
an important parameter via this mechanism. 

We emphasize an incompleteness of the present survey to- 
wards higher densities, at which electron captures in the ashes 
become important. These shift the burning products to neutron- 
rich isotopes, favoring 58 Ni instead of 56 Ni. This effect would 
be taken into account in our postprocessing procedure, how- 
ever electron captures may also become dynamically important 
with increasing central densities, since they reduce the electron 
pressure in the ashes. Unfortunately, this effect could not con- 
sistently be modeled in the current study. The explosion model 
assumes Y e = 0.5. Effects of higher central densities will be 
addressed in forthcoming investigations. 

7.1. Comparison with one-dimensional models 

The effect of a variation in the carbon mass fracti on of the pro- 
genitor on the produced 56 Ni mass was studied bv lHoflichetai] 
(1998). They applied a one-dimensional delayed detonation 
model. For a central density of 2.6 x 10 9 gcirT 3 , solar metal- 
licity and a presumed deflagration-to-detonation transition at a 
density of 2.7 x 10 7 g crrT 3 they calculated a model with a C/O 
ratio of 1/1 (DD21c in their notation) and a model with C/O re- 
duced to 2/3. Here they find a decrease of the produced M( 56 Ni) 
of about 14%. Assuming that a transition to detonation at such 
low densities as applied here does not alter the production of 
iron peak elements, a comparison with our models is possible. 
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However, the results o f lH6flichetal] Jl998) are in contrast with 
ours. This may be mainly due to the fact that the modeling of 
the correct implications of the C/O ratio on the explosion re- 
sults requires an accurate description of the multidimensional 
effects that dominate the flame propagation. 

Brav o et al] (^993) investigated the impact of the igni- 
tion density on the 56 Ni production for one-dimensional de- 
flagration models. For models with a central density of 2.5 x 
10 9 g cirT 3 (R2 in their notation) and 4.0 x 10 9 g cirT 3 (R4) they 
find differences of about 7% which is in good agreement with 
our results. 

Although our results regarding the change in 56 Ni produc- 
tion varying the metalli city are in good agr eement with the an- 
alyti cal prediction by Timm es et alJ (12003) and with the study 
bv iTrava glio et al] d2005l) . thev disagree with the findings of 
iHoflich et all (1998). They report an only ~5% effect on the 
56 Ni production changing the metalli city from 0.1Z o to 1OZ . 
Contrary to this, the result of llwamoto et alJ lfl999) that an in- 
crease of the metallicity from zero to solar decreases the 56 Ni 
production for about 8% is consistent with our results. 



7.2. Cosmological significance 

In order to discuss the cosmological significance of our results, 
we take on a very simplistic view on the mecha nism of the light 
curve. Following "Arnett's rule" ( lArnettLfl982h we assume that 
the total mass of 56 Ni immediately determines the peak lu- 
minosity of the SN la event. Furthermore we assume that a 
larger energy released in the explosi on leads to a m ore rapid 
decline of the light curve jPinto & Eastman! 12000). It has to 
be noted, however, that this may only be a second-order effect 
in the light-curve shape. The main parameter h ere is the opac- 
ity gi ven by the distribution of heavy elements ( Ma zzali et all 
l200ll) . This can only be adequately addressed in detailed syn- 
thetic light curve calculations and will be ignored here. 

In the context of this simplification, we note that none of the 
tested parameters reproduces the peak luminosity-light curve 
shape relation by lowering the produced 56 Ni mass accom- 
panied by an increased energy release. While the carbon-to- 
oxygen ratio of the progenitor has little effect on the peak lu- 
minosity, it could alter the width of the light curve. The oppo- 
site holds for the progenitor's metallicity. Here, the peak lumi- 
nosity can be vary by about 20%, but the explosion dynamics 
is unaffected. The central density prior to the ignition changes 
both the 56 Ni production and the energy release. Unfortunately 
our study is incomplete here. At higher values for the central 
density, the produced 56 Ni mass could decrease due to electron 
captures while the energy release may still increase. This has 
to be tested in forthcoming studies. 

Another aspect is that we have ignored the interrela- 
tion of the parameters by stellar evolution here. This, how- 
ever, predicts a lowe r C/O ratio for higher metallicities (cf. 
Umed aetali [l999b). The effects of both parameters in this 
combination may possibly reproduce the trend of the peak 
luminosity-light curve shape relation. 

The final conclusions on the cosmological significance of 
the variations in the explosions found in the present study need 



to be drawn on the basis of synthetic light curves derived from 
our models. This is subject of a subsequent publication. 
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Table 6. Synthesized mass (M Q ) for different models as plotted in Figs. II lllT2l and !13l 
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-04 


34 S 


1 


.0709E- 


■03 


1 


.1193E- 


■03 


1 


.2621E- 


■03 


1 


.3314E- 


■03 


5 


4058E- 


-04 


3 


.0902E- 


-03 


36 s 


5 


.0047E- 


-07 


2 


.4312E- 


■07 


1 


.5159E- 


-07 


2 


.8158E- 


■07 


3 


.6821E- 


-08 


4 


.711 1E- 


-06 


35 CI 


5 


. 042 IE- 


■05 


5 


.0890E- 


■05 


4 


.9569E- 


■05 


5 


.2873E- 


■05 


3 


.0618E- 


-05 


6 


.0507E- 


-05 


37 C1 


1 


.3557E- 


■05 


1 


.023 IE- 


■05 


7 


.5430E- 


06 


1 


.1730E- 


■05 


7 


.0770E- 


-06 


1 


.5845E- 


-05 


36 Ar 


4 


.7197E- 


■03 


3 


.5313E- 


■03 


2 


.6786E- 


■03 


3 


.9809E- 


■03 


3 


.5923E- 


-03 


2 


.6380E- 


-03 


3S Ar 


5 


.9756E- 


-04 


4 


.8536E- 


■04 


3 


.5572E- 


■04 


5 


.7181E- 


-04 


2 


.2226E- 


-04 


1 


.5216E- 


-03 


4.0 a 

4U Ar 


6 


.3966E- 


09 


2 


.3371E- 


09 


1 


.291 IE- 


09 


2 


.6198E- 


09 


2 


.6172E- 


-10 


5 


.5970E- 


-08 




4 


.1443E- 


■05 


2 


.4448E- 


■05 


1 


.1254E- 


■05 


2 


.8140E- 


■05 


1 


4635E- 


-05 


4 


.1644E- 


-05 


40 K 


1 


.5893E- 


-08 


1 


.1857E- 


-08 


8 


.8932E- 


-09 


1 


.2172E- 


-08 


3 


.5607E- 


-09 


1 


.5243E- 


-08 


41 K 


3 


.3541E- 


06 


1 


.6253E- 


06 


6 


.437 IE- 


■07 


1 


.8679E- 


06 


1 


.0319E- 


-06 


2 


.3920E- 


-06 


A ft „, 

40 Ca 


4 


.1245E- 


■03 


2 


.9254E- 


■03 


2 


.3081E- 


■03 


3 


.2055E- 


■03 


3 


.0169E- 


-03 


2 


.2343E- 


-03 


42 Ca 


2 


.0467E- 


■05 


1 


.1405E- 


■05 


5 


.0761E- 


06 


1 


.3292E- 


■05 


4 


.671 IE- 


-06 


3 


.7149E- 


-05 


43 Ca 


3 


.7921E- 


-08 


3 


.3294E- 


08 


3 


.7149E- 


08 


3 


.6033E- 


-08 


3 


.7361E- 


-08 


5 


.5808E- 


-08 


44 Ca 


2 


.7296E- 


06 


2 


.8433E- 


06 


3 


.2414E- 


-06 


2 


.9370E- 


06 


3 


.3912E- 


-06 


1 


.3942E- 


06 


46 Ca 


1 


.1429E- 


11 


3 


.9578E- 


-12 


1 


.7961E- 


-12 


4 


.6778E- 


■12 


3 


.6129E- 


-13 


6 


.5306E- 


11 


4S Ca 


5 


.1102E- 


16 


4 


.0742E- 


■17 


8 


.4593E- 


-18 


4 


.723 IE- 


-17 


7 


.0210E- 


19 


3 


.8800E- 


-14 


45 Sc 


9 


.5903E- 


-08 


4 


.1321E- 


08 


1 


.8759E- 


-08 


4 


.5593E- 


-08 


2 


.8767E- 


-08 


5 


.8257E- 


-08 


46 T j 


8 


.5527E- 


06 


4 


.8903E- 


06 


2 


.2655E- 


-06 


5 


.609 IE- 


06 


2 


.1275E- 


-06 


1 


.3952E- 


-05 


4/ Ti 


2 


.0932E- 


-07 


2 


.0933E- 


-07 


2 


.1956E- 


■07 


2 


.1773E- 


-07 


2 


.0015E- 


-07 


4 


.9531E- 


-07 


48 Tl 


1 


.5686E- 


-04 


1 


.4112E- 


-04 


1 


.3379E- 


-04 


1 


.3985E- 


-04 


1 


.5444E- 


-04 


1 


.0256E- 


-04 


4Q m . 

4J Tl 


6 


.2734E- 


06 


5 


.3616E- 


06 


4 


.8099E- 


06 


5 


.1362E- 


06 


4 


.1781E- 


-06 


7 


.1369E- 


-06 


5() Ti 


4 


.0498E- 


-10 


3 


.9554E- 


-10 


4 


.0308E- 


-10 


3 


.3934E- 


1 1 


3 


.5590E- 


-10 


1 


.1353E- 


-09 


50y 


1 


.6127E- 


09 


1 


.5420E- 


09 


1 


.8142E- 


09 


4 


.2838E- 


■10 


1 


.240 IE- 


-09 


6 


.9824E- 


-09 


51y 


2 


.1755E- 


■05 


1 


.8974E- 


■05 


1 


.7403E- 


■05 


1 


.3477E- 


■05 


1 


.5028E- 


-05 


3 


.3895E- 


-05 


50 Cr 


1 


.2996E- 


-04 


1 


.2398E- 


-04 


1 


.0947E- 


■04 


9 


.1295E- 


■05 


8 


.6254E- 


-05 


2 


.6088E- 


-04 


^9 ^ 

52 Cr 


2 


.7097E- 


-03 


2 


.4212E- 


03 


2 


.2432E- 


03 


1 


.7785E- 


■03 


2 


.5029E- 


-03 


2 


.4744E- 


-03 


Cr 


c 

J 


.Zoo ih- 


-U4 


1 

4 


.ol /9b- 


-U4 


4 


,433ob- 


-04 


2 


.4196b- 


-U4 


4 


.2566b- 


-U4 


c 
6 


-7 A n cc 

. /4/jb- 


-U4 


54 Cr 


2 


.8910E- 


-08 


2 


.9013E- 


■08 


2 


.8687E- 


-08 


5 


.3360E- 


-11 


2 


.785 IE- 


-08 


3 


.6635E- 


-08 


55 Mn 


6 


.0348E- 


■03 


5 


.6617E- 


■03 


5 


.2371E- 


-03 


2 


.8042E- 


-03 


5 


.1364E- 


-03 


7 


.7788E- 


-03 


54 Fe 


6 


.8344E- 


-02 


6 


.4350E- 


-02 


5 


.9272E- 


-02 


2 


.9665E- 


-02 


5 


.9137E- 


-02 


8 


.7595E- 


-02 


56p e 


3 


.2193E- 


-01 


3 


.2643E- 


■01 


3 


.281 IE- 


■01 


2 


.8657E- 


01 


3 


.3773E- 


-01 


2 


.8308E- 


-01 


57 Fe 


1 


.2662E- 


■02 


1 


.3360E- 


■02 


1 


.3723E- 


■02 


1 


.0757E- 


■02 


1 


.2245E- 


-02 


1 


.5385E- 


-02 
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Table 6. continued. 



Species model 







2_7_2 






2_2_2 






2_3_2 






7_2_2 






2_2_1 






2_2_3 




S8 Fe 


6 


4899E- 


■06 


6. 


6366E- 


■06 


6 


.7433E- 


■06 


1 


.2656E- 


■07 


6. 


.5205E- 


■06 


7. 


.2146E- 


06 


59 Co 


5. 


8039E- 


■04 


6. 


.1402E- 


■04 


6 


4570E- 


■04 


2 


.0260E- 


■04 


6. 


2006E- 


■04 


6. 


1134E- 


04 


58 Ni 


6 


.6717E- 


■02 


1. 


0699E- 


•02 


7. 


.3404E- 


•02 


4 


.5109E- 


■02 


6. 


4330E- 


■02 


9 


.2729E- 


02 


60 Ni 


6 


.3865E- 


■03 


6. 


9060E- 


■03 


7 


.3461E- 


■03 


1 


.9357E- 


■03 


1. 


3971E- 


■03 


6. 


0938E- 


03 


61 Ni 


6. 


3095E- 


■05 


8. 


1359E- 


•05 


9 


.5324E- 


■05 


7 


.6708E- 


■05 


8. 


.5555E- 


■05 


4. 


3643E- 


05 


62 Ni 


5. 


.2705E- 


■04 


6. 


9753E- 


■04 


8 


.2941E- 


■04 


6 


.6740E- 


04 


5. 


8056E- 


■04 


6. 


.1748E- 


04 


64 Ni 


3. 


4339E- 


■11 


3, 


,8872E- 


•11 


4 


.2017E- 


■11 


5 


.6274E- 


■14 


3. 


,7882E- 


■11 


4. 


3591E- 


-11 


63 Cu 


7. 


8178E- 


■07 


1 


.0014E- 


■06 


1 


.1137E- 


06 


7 


.5785E- 


■07 


1 


4388E- 


■06 


7. 


3209E- 


07 


65 Cu 


1 


.8541E- 


■07 


2. 


.504 1 E- 


■07 


2 


.9743E- 


■07 


2 


.4756E- 


■07 


2. 


8771E- 


■07 


9 


1714E- 


08 


64 Zn 


4 


1170E- 


06 


5. 


403 IE- 


■06 


6 


.3524E- 


■06 


5 


.5728E- 


06 


1. 


0068E- 


■05 


1. 


.2144E- 


06 


66 Zn 


5, 


.3276E- 


■06 


7. 


.1634E- 


■06 


8 


.5197E- 


06 


6 


.8650E- 


■06 


6. 


0776E- 


06 


5. 


.0221E- 


06 


67 Zn 


2. 


.1888E- 


■12 


2. 


.7503E- 


■12 


3 


.2022E- 


■12 


1 


.0714E- 


12 


2. 


.5332E- 


■12 


2. 


.9738E- 


12 


68 Zn 


2. 


4497E- 


■09 


3, 


4762E- 


09 


4 


.2682E- 


■09 


2 


.9172E- 


09 


2. 


4247E- 


09 


6. 


.0432E- 


-09 


™Zn 


L 


6091E- 


■19 


2. 


1266E- 


■19 


2 


.3458E- 


■19 


1 


.4498E- 


■23 


2. 


0416E- 


■19 


2. 


5513E- 


-19 



F. K. Ropke et al.: Type la supernova diversity in three-dimensional models, Online Material p 4 
Table 7. Synthesized mass (M Q ) of radioactive species in different models 



Species model 





2_1_2 




2_2_2 




2_3_2 






1_2_2 




2_2_7 






2_2_3 




22 Na 


5.2892E 


08 


9.2173E 


08 


1.0290E- 


07 


8 


0645E 


08 


9.9089E 


08 


4 


1895E 


08 


26 Al 


1.3078E 


06 


1.1027E- 


06 


8.0954E- 


07 


1 


0024E 


06 


4.71 10E- 


07 


2 


8228E 


06 


36 C1 


7.6022E 


07 


7.2749E 


07 


6.8645E- 


07 


8 


1239E 


07 


2.8832E 


07 


1 


1296E 


06 


39 Ar 


4.3326E 


09 


2.653 1E- 


09 


1.7073E- 


09 


2 


5865E 


09 


5.3153E 


10 


1 


3045E- 


08 


40 K 


1.5893E 


08 


1.1857E- 


08 


8.8932E- 


09 


1 


2172E- 


08 


3.5607E 


09 


1 


5243E- 


08 


41 Ca 


3.3488E 


06 


1.6224E 


06 


6.4171E- 


07 


1 


8647E 


06 


1.031 1E- 


06 


2 


3714E- 


06 


44 Ti 


2.7207E 


06 


2.8370E- 


06 


3.2374E- 


06 


2 


929 IE 


06 


3.3876E 


06 


1 


3530E 


06 


48y 


2.7174E- 


08 


1.6707E 


08 


1.0033E- 


08 


1 


9273E- 


08 


1.0073E- 


08 


4 


0923E- 


08 


49y 


6.0178E- 


08 


5.0885E 


08 


5.9430E- 


08 


3 


3875E 


08 


3.0780E 


08 


2 


8737E- 


07 


53 Mn 


2.0549E 


04 


1.9908E 


04 


1.9552E- 


04 


1 


6387E 


05 


1.9147E- 


04 


2 


7557E- 


04 


55p e 


2.0833E 


03 


2.0297E 


03 


1.9799E- 


03 


2 


1947E 


04 


1.9797E 


03 


2 


3335E- 


03 


59 Fe 


8.7600E 


13 


9.1399E 


13 


9.483 1E- 


13 


8 


2378E 


17 


8.7169E 


13 


1 


1308E- 


12 


60p e 


4.7700E 


15 


5.2457E 


15 


5.5233E- 


15 


1 


2764E 


19 


4.9474E 


15 


6 


8064E 


15 


56 Co 


1.0058E 


04 


9.8004E 


05 


9.3479E- 


05 


4 


5078E 


05 


9.2616E 


05 


1 


248 IE- 


04 


57 Co 


1.077 IE- 


03 


1.0724E- 


03 


1.0677E- 


03 


1 


1285E 


04 


1.053 IE 


03 


1 


1640E- 


03 


60 Co 


2.5561E 


10 


2.7429E 


10 


2.8596E- 


10 


1 


8310E 


13 


2.6598E 


10 


3 


1485E- 


10 


56 Ni 


2.9992E 


01 


3. 046 IE 


01 


3.0648E- 


01 


2 


8599E 


01 


3.1625E 


01 


2 


5917E- 


01 


57 Ni 


1.1578E 


02 


1.2281E- 


02 


1.2648E- 


02 


1 


0643E 


02 


1.1183E- 


02 


1 


4214E- 


02 


59 Ni 


4.6795E 


04 


4.78 14E- 


04 


4.891 1E- 


04 


8 


0744E 


05 


4.6775E- 


04 


5 


2852E- 


04 


60 Ni 


5.1163E 


03 


5.2800E 


03 


5.4430E- 


03 


2 


3015E- 


04 


5.2058E 


03 


5 


6263E 


03 


61 Ni 


2.8036E 


06 


2.9852E 


06 


3.1293E- 


06 


2 


4920E 


08 


2.9360E 


06 


3 


2136E 


06 


62 Ni 


5.3714E 


06 


5.7961E- 


06 


6.1280E- 


06 


8 


1844E 


09 


5.6628E 


06 


6 


4208E 


06 


63 Ni 


2.2585E 


11 


2.5002E 


11 


2.6783E 


11 


3 


4056E 


15 


2.4055E 


11 


2 


9699E 


11 



